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Abstract 

Cilia are motile biological appendages that are driven to bend by internal shear stresses between 
tubulin filaments. A continuum model of ciliary material is constructed that incorporates the 
essential ciliary constraints: (1) one-dimensional inextensibility of filaments, (2) three-dimensional 
incompressibility, and (3) shear strain only along filaments. This hypothetical ciliary material 
combines one- and three-dimensional properties in a way that makes it a natural and flexible model 
for how real cilia convert nanoscopic shear stress into motility on a much larger scale. Without 
reference to the evolving shape of the cilium, conventional continuum mechanics applied to this 
hypothetical material leads to the standard model of ciliary dynamics, but with one additional 
term, required by constraints (2) and (3) above, a model-independent coupling of shear and twist 
in general ciliary motion. 
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I. INTRODUCTION 



The beating of cilia, or eukaryotic flagella, has been the subject of modelling studies for 
over fifty years. The hydrodynamics of propulsion is one aspect of these studies, concerned 
with Stokes flows past the cilium at a mesoscale of tens of microns [l2]. This scale 

is set by the length of the cilium. The mechanism of motility is another aspect, concerned 
with the dynein mediated sliding of microtubule doublets on each other at a nanoscale of 
tens of nanometers f| ^Qjfs] 10]. This scale is set by the diameter of the cilium or even 
by the size of the motor proteins. It is the aim of this paper to bridge the gap between these 
scales. 

The remarkable discrepancy between these length scales poses a problem for modelling. 
To put it simply, at the mesoscale the cilium looks one-dimensional, but at the nanoscale 
it looks three-dimensional. The discrepancy is most telling when the ciliary beat itself is 
three-dimensional. In modelling the cilium at the mesoscale, the emphasis is on the exterior 
of the cilium regarded as a space curve. The terms of the description are dictated by this 
curve, whether in the Frenet description of Gueron and Liron's three-dimensional ciliary beat 
model [ll], or in the 'body co-ordinates' of a later version by Gueron and Levit-Gurevich 
[l^ . The three-dimensional nanoscopic models of Hines and Blum [3] Q] [2 ] ID], on the other 
hand, use body coordinates truly anchored in the three-dimensional cilium, and moving with 
it, but these authors consider only static stress-strain relationships and do not address the 
three-dimensional mesoscale dynamics of the cilium. Ultimately one would like to bridge 
the gap between these scales so that observations of behavior at the mesoscale, which are 
increasingly accessible through high speed video microscopy, could be simply related to 
geometry at the nanoscale. 

This article is essentially an abstraction from the sliding filament models of Hines 
and Blum. Although these authors work on detailed three-dimensional structures at the 
nanoscale, they consistently summarize the results in effective parameters for mesoscale 
models. In this article I determine what must be true of all such models, given only the 
constraints of sliding filament models in general. It turns out that the constraints at the 
nanoscale show up in dynamics at the mesoscale, and this is a model-independent result. 

The plan of the paper is as follows. Section [Til gives an abstract characterization of a three 
dimensional material that I call 'ciliary material,' subject to the constraints of the sliding 
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filament model: it is incompressible, inextensible along one-dimensional "filaments," and 
admits shear strain only longitudinally along the filaments. This is an attempt to capture, 
in a continuum model, the properties of the cilium. These are the constraints consistently 
invoked by Hines and Blum, but without reference to nanoscopic details. These constraints 
tightly couple longitudinal shear and bending strain, the basic mechanism of ciliary action. 
Such ciliary material might also model other biological structures, like auditory hair cells, 
where this coupling has a different function. 

Section HTT1 determines how the ciliary material can move, subject to its constraints, what 
I call ciliary flow. Remarkably, although it is a three-dimensional material, its motion is 
parameterized by the motion of any single one- dimensional filament, and this one filament 
can move arbitrarily. Thus the hypothetical ciliary material, motivated by the structure of 
cilia, combines, in an unexpected way, three-dimensional and one-dimensional properties in 
a manner that is reminiscent of real cilia. 

The transition in Section IIVI from the nanoscopic properties of three-dimensional ciliary 
material to the mesoscale dynamics of an essentially one-dimensional cilium is a very short 
step. The treatment there looks one dimensional, but the description retains information 
about the three dimensional nanostructure because of the peculiar properties of ciliary ma- 
terial. The most significant result is a model-independent coupling between shear and twist 
in three dimensional ciliary beats, a phenomenon discussed by Hines and Blum in specific 
models, but not recognized by them as a generally necessary feature of every sliding filament 
model. 

The geometrical methods in Sections [TT] and III II are not used very often in physics, but 
they are indispensable here, and I hope any reader who is unfamiliar with them will find 
that this application is a useful introduction to them in a relevant context. The books by 
Schutz [13I and Frankel [ljj] are good references. 

II. CILIARY MATERIAL 

In terms of smooth coordinates (a; 1 , a; 2 , x 3 ) in space one can describe the deformation of 
any material by the trajectories of its constituent particles, solutions of equations of motion 




(1) 
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where 

V = V% (2) 
is the velocity vector field, and t is time. Metric relations among particles are given by 

ds 2 = gijdx l dx^ (3) 

where gij is the Euclidean metric. Coordinatize the material object by Lagrangian coordi- 
nates convected by the flow, i.e., let every material point keep the same coordinates that 
it had originally. The changing metric relationship of material points is then expressed en- 
tirely by the change in the metric components g^, and the rate of change is given by the Lie 
derivative 

99jk ~ Vg jk + g([d jt V],d k ) + gidj, [d k , V}) (4) 



dt 

Here [ , ] is the Lie bracket of vector fields. 

As an example of such a computation, let V be the shearing flow 

V = ySd x (5) 

in the plane, where S is a constant. Then if x and y are initially Cartesian coordinates, 
the metric tensor in convected coordinates at time t is given by the Lie- Taylor series (which 
terminates in this case) 

t 2 ( 1 tS \ 

g + t£ v g + —£ v £ v g = , (6) 

21 \tS l + t 2 S 2 ) 

and evaluating at t = 1 one has the Euclidean metric 

(i s \ 

in skew coordinates corresponding to constant shear strain S. 

The flow V of Eq. (J5]) is an example of a ciliary flow, and the form of the metric tensor 
g in Eq. (JTj) encodes the result. The filaments, parallel to the x-axis and labeled by the y 
coordinate, slide on each other inextensibly, as one sees in the metric component g\\ = 1. 
The resulting shear strain is visible in the off-diagonal component gu = S. Finally, det g = 1 
means that the flow is incompressible. The three dimensional generalization of this form to 
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a general ciliary configuration is 

(lS T \ 

S 1 + S 2 ST 



(8) 



\T ST 1 + T 2 J 

where now S and T, two independent components of shear strain, depend on Lagrangian 
coordinates (x, y, z). Here x is arclength along filaments, up to an additive constant reflecting 
the arbitrariness of choosing an origin in each filament, and the coordinates (y, z) label the 
filaments. It should be noted that this is not the most general outcome of an incompressible 
flow of the filaments, because the only motion that has been allowed to them is sliding 
longitudinally along each other (and of course bending). These are exactly the constraints 
of the sliding filament models of Hines and Blum, encoded in the form of g. 
Associated with the metric tensor g is an orthonormal frame field 

ei = d x (9) 
e 2 = d y - Sd x (10) 
e 3 = d z -Td x . (11) 

where e\ is everywhere tangent to filaments. This is not, however, a coordinate frame, 
because the vector fields do not commute as differential operators, in general. 

The metric g determines the shape of each filament up to a rigid motion, because for 
each (y, z) it determines the filament's curvature k and torsion r as a function of x (see 
Appendix [A]). These are the Frenet data for the filament as a space curve. One has 



K=jS* + I* (12) 
T <? — <7 T 

J -xx LJ x LJ xx J -x /10\ 

r = (13) 

K 

where the subscript x denotes the derivative d x along the filament. The Frenet equations for 
the filaments then determine the space curves R(x) and an orthonormal frame {T, N, B} at 
each point of each curve according to 

R x = f (14) 

f x = kN (15) 

N x = - k T + tB (16) 

B x = -tN (17) 



In the context of ciliary geometry it is more natural to describe the filaments in terms of 
the shear strains S and T directly than to translate them into the Frenet language. This 
is an alternative, equivalent description of space curves. It is even possible to translate the 
other way, from the Frenet description to the ciliary description, inverting the relations in 
Eqs. (tH-dH, 

S x = k cos (f) (18) 
T x = Ksin0 (19) 

0* = T (20) 

In this formulation the space curves R(x) are the solutions of 

R x = ei (21) 

d x e\ = Kcos0e2 + Ksin0e3 (22) 

d x e2 = — /tcos0ei (23) 

d x es = —Ksmcfiei (24) 

0* = r (25) 

There is an arbitrary constant in the angle since, unlike N and B in the Frenet picture, 
e 2 and e3 are only determined up to a global rotation about e±. This amounts to a gauge 
freedom in the ciliary description. An advantage of the ciliary description is that it remains 
well defined where the curvature vanishes, a nuisance in the Frenet description. 

The shear strains S and T cannot be arbitrary functions, because g is the Euclidean 
metric, even if it is expressed in peculiar coordinates. Its associated Riemannian curvature 
tensor must therefore vanish identically, and hence S and T must obey the following identities 
(see Appendix [B]) 

= d x (e 2 S) = d x (e 3 S) = d x (e 2 T) = d x (e 3 T) (26) 
= e 3 S - e 2 T . (27) 

The second of these equations can be recognized as 

[e 2 ,e 3 ]=0. (28) 

This says that e 2 and e 3 together form an integrable distribution, and thus there exist 
surfaces normal to the filaments, at least locally. Moreover, since e 2 and e 3 are now an 



orthonormal coordinate frame field on these surfaces, they are Euclidean planes, which I 
will call normal planes. 

By allowing the filaments to reptate along their length inextensibly, one can bring Eq. (126]) 
to the simpler form 

= e 2 S = e 3 S = e 2 T = e 3 T (29) 

that is, the strains S and T can be made constant in normal planes. Thus the configuration 
of filaments in a neighborhood of a given filament is entirely determined by that filament. 
The given filament determines the direction of neighboring filaments, since it determines 
the normal planes, and it determines their curvatures and torsions, since it determines the 
values of S and T. The geometrical meaning of Eqs. (p26]) - (l2"Tj) is that these two potentially 
conflicting descriptions are consistent. Thus the ciliary geometry is almost rigid. The 
freedom that it possesses corresponds essentially to a single free space curve. 

I turn now to the motions possible in ciliary matter, ciliary flows. Since one of those 
motions is reptation of filaments, I will also justify the assertions of the preceding paragraph. 



III. CILIARY FLOWS 



A ciliary flow, by definition, must maintain the form of Eq. ([8]) even as g evolves according 
to Eq. (jl]). If V is the flow with components (a, (3, 7), namely 

V = ad x + (3d y + 1 d z , (30) 

then for V to be ciliary (a,/?, 7) must obey the conditions (see Appendix ICl) 

= a x + Sf3 x + T lx (31) 

= e 2 (3 (32) 

= e 37 (33) 

= e 2l + e 3 (3 (34) 

= {e 2 r r)x - PxT x + "i x S x (35) 

= e 2 e 2 7 (36) 

= e 3 e 3 (3 (37) 

where the partial derivative d x is indicated by the subscript x. Eqs. (I36]) - (I37]) are not 

independent of the others, since for example e 3 e 3 /3 = — e 3 e 2 7 = — e 2 e 3 7 = 0. In fact all 



higher derivatives of (3 and 7 in the normal planes vanish by this argument, so that (3 and 7 
restricted to a normal plane can only be affine linear functions. I return to this consideration 
below. 

The simplest nontrivial ciliary flow is reptation, a x = (3 = 7 = 0. By Eq. fll]) the shear 
strains change under this flow at the rate 

S t = (x y + aS x (38) 
T t = a z + aT x , (39) 

where subscripts indicate partial derivatives with respect to the corresponding variable. It is 
possible to construct a reptation a that alters S and T in such a way that Eq. (129]) holds in 
one normal plane, at least locally (see Appendix [D]). Then by Eq. fl26|) S and T are constant 
in every normal plane along the filaments. This proves the assertions made at the end of 
the last section. I will now assume that S and T obey the simpler Eq. ( 1291) . since this can 
always be arranged by a reptation. 

The ciliary conditions require f3 and 7 to be affine linear functions in the normal planes. 
To see that there exist non-trivial solutions to these conditions, imagine specifying (3{x) and 
7(x) arbitrarily along one filament. Integrating Eq. (|35|) determines e27 = —e^P along the 
filament up to an arbitrary global constant. These are all the data required to extend j3 
and 7 in each normal plane as affine linear functions. This extension continues to satisfy 
Eq. ( 135]) on neighboring filaments because of the easily verified identities 

e 2 [(e 2 7)x - T X (3 X + S x j x ] = 2S x [(e 2 ^) x - T X (3 X + S x ^ x ] (40) 
^{{e 2 i) x - T X (3 X + S x j x ] = 2T x [(e 2 j) x - T X (3 X + S x j x ) (41) 

Repeated differentiation of Eq. ( 135]) using e 2 and shows that if (e 2 j) x — T X (3 X + 
vanishes on a filament, then all its normal derivatives also vanish there. Thus given only 
that it is represented by its Taylor series, it is constant, and hence zero, and the constructed 
solution obeys all the ciliary conditions. V is a non-trivial ciliary flow, determined by its 
values on one filament. 

The conditions Eqs. (|3T|)-(l37j) confirm that the ciliary configuration is entirely determined 
by one filament. If one tries to move the filaments, one can specify f3(x) and 7(2) on only 
one given filament. The ciliary conditions then determine a(x) along that filament up to a 
constant (a reptation). They further determine (3 and 7 almost uniquely as linear functions 
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in normal planes. Neighboring filaments intersect the normal planes (each plane labelled by 
x, its intersection with the given filament), and in the course of the motion these intersection 
points rotate about the given filament at an angular velocity 

u){x) = (e 2 l)(x) (42) 

which is not arbitrary but is determined up to a global constant by (3 and 7 according to 
Eq. (1351) . (The undetermined constant in u describes a global twisting rotation in which every 
normal plane rotates at the same rate, a motion which is possible but probably not physically 
relevant.) The remaining freedom in a corresponds to reptations of neighboring filaments, 
also probably not physically relevant. Thus, to summarize, any filament determines the 
configurations of its neighbors, and any motion of that filament determines the motion of 
its neighbors. It is worth noting what was not obvious a priori, that non-trivial motions of 
the filaments are possible, that is, the ciliary material is not completely rigid. One filament 
can move arbitrarily (but inextensibly) , and all the others must follow it. 



IV. CILIARY DYNAMICS 



The surprisingly one-dimensional character of ciliary matter means that its restriction 
to modelling the cilium is immediate. Represent the cilium as a thin cylinder of ciliary 
matter, with length L and radius p, all dynamical quantities now depending just on x, the 
coordinate along one given filament. Let the cilium have elastic moduli, so that the energy 
of a configuration is given by the shear energy and bending energy 

E = \ J\{S - Ff + (T — G) 2 ] dx+^ J\sl + T 2 ) dx , (43) 

where [i is the shear modulus and k c is the bending modulus (with dimensions appropriate 
to one dimension, not three). Here F and G are target shear strains, such that the shear 
energy would be minimal if the shear strains S and T could relax to F and G. The operation 
of the ciliary "engine" would be to make F and G functions of time, so that the equilibrium 
becomes a moving target. Asymmetries in the nanoscopic structure of the cilium could be 
built into this expression, which is taken symmetrical here for illustration. 
Conjugate to shear strains S and T are internal stresses 

$ = 5E/5S = fi(S — F) — k c S xx (44) 
= SE/6T = /i(T — G) — k c T xx (45) 
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which, by the above mentioned one-dimensionality, have the units of force. 

Continuing this variational approach to the dynamics, one finds the generalized forces 
conjugate to the displacement of a cilium along V = ad x + j3d y + 7<9 2 of Eq. (I3"U|) . together 
with the required rotation uj of Eq. (T421 in normal planes. In the orthonormal basis, V takes 
the form 

V = (a + (3S + 7 T)ei + (3e 2 + 7 e 3 , (46) 

a representation that becomes increasingly appropriate as one moves up to the mesoscopic 
scale. One notes that (3 and 7 are the components of velocity normal to the cilium, but that 
they also contribute to the tangential velocity if the shear strains are non-zero. This is a 
residual piece of three dimensional information in the one-dimensional description. 

From Eq. (j4j), imposing Eqs. (l3TT) - (l37l) . the strains under the flow V change at the rate 



S t = (a + pS + ir)S x + P X + e 2 a + coT (47) 

T t = (a + (3S + jT)T x + lx + e 3 a - uS (48) 

The terms €20. and represent reptations within the cilium, a possibility that is ignored 
from now on. The constraints 

= a x + S(3 X + T lx (49) 

= u x - T x f3 x + S x "i x (50) 

are handled with Lagrange multipliers A and v in the expression 

E' = E + [ \(a x + Sp x + Tj x ) dx + [ v(lo x — T X (3 X + S x ^ x ) dx (51) 
Jo Jo 

Then the generalized interior forces on the cilium are 

F a = -50 /5a = -®S X - VT X + X x (52) 

Fp = -5E'/5(3 = -§SS X - t>ST x + $ x . + (XS) X - {vT x ) x (53) 

F 7 = -5E'/5 1 = -$TS X - VTT X + V x + (XT) X + (uS x ) x (54) 

= -5E'/5u = v x - $T + (55) 



Generalized external forces on the cilium, due to the fluid medium of viscosity r] in which 
it is immersed, can be found from the dissipation function D for the corresponding Stokes 
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flow. A thorough analysis of this problem has been done by Gueron and Liron [4]. For 
simplicity I take just the leading term indicated by their analysis, 



D 



o 



— {a + pS + ^T) +-2"(/3 +7) + -2~w 



c?x (56) 



where the C's are constant resistance coefficients. The viscous force conjugate to displace- 
ment by a is then —5D/5a, etc. Requiring that elastic forces and viscous forces balance, 

5a 5a ^ ^ 
etc., leads to the surprisingly simple dynamical laws connecting the motion of ciliary matter 
to its internal stresses, 

C T (a + (3S + 7 T) = -<S>S X - *T X + \ x (58) 

C N (3 = + \S X - {vT x ) x (59) 

C Nl = V x + XT X + (uS x ) x (60) 

Cjjj = v x -$T + tyS (61) 

Under this flow the strains change according to Eqs. (j47l) - fj48l . The Lagrange multipliers 
can now be determined from the constraints Eqs. (|49|) - (|50|) . Substituting the solutions in 
Eqs. f l58|) - fl6T]) gives the equations for A and u, 

Kx = + ($S X + *T X ) X + ^{S x ^ x + T x * x ) + ^tk 2 v (62) 

c 

v xx = ($T - $?S) X - — - S x ty xx - T x $ xx + Xtk 2 + T x (vT x ) xx + S x (uS x ) xx , (63) 

where I have used the expressions for curvature k and torsion r first derived in terms of S and 
T in Eqs. (JT2l)-(fT3l. In the limit as C^/Cn — > 0, the equation for v simplifies considerably, 
with solution 

v= / 2 ($T - VS) dx + c x x + c 2 (64) 
Jo 

The solution for u is already explicit in terms of (3 and 7 from the constraint Eq. fl50|) . 

oj = (T X (3 X - S x >y x ) dx + const. (65) 
Jo 

With one important exception, the above dynamical laws for ciliary matter are precisely 
the usual phenomenological laws for three-dimensional motions of cilia, as derived by Gueron 
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and Liron . It is remarkable that they emerge here without any appeal or reference to the 
shape of the cilium (apart from the hydro dynamic interaction in a more accurate dissipation 
function, Eq. (1351) ). but purely as a consequence of nanoscopic processes. The occurrence 
of k and r above is just an abbreviation for certain combinations of derivatives of S and T 
that appeared. Of course the large scale shape of the cilium could be reconstructed at any 
time from Eqs. (fT8])-([25l). 

The dynamical laws have been expressed here in terms of the force vector 

F = Aei + $e 2 + ^e 3 (66) 



11] were expressed in terms of the Frenet basis T, N, B, where 



but the vector quantities in 
F took the form 

F = -F T f - F N N - F B B . (67) 
From Eqs. (|T8l)-(l22l) one has the transformation that connects these two descriptions, 

f = ei (68) 
N = cos e 2 + sin e3 (69) 
B = — sin0e2 + cos0e3 (70) 



with <p x = t, the torsion of the curve. Using also Eqs. ( 112l) - (lT3l) . it is straightforward 
to verify that the dynamical laws of ciliary matter, as derived here, agree with standard 
phenomenology in every respect but one: in ciliary matter there is an additional constraint, 
and a corresponding additional term in the laws. 

The additional constraint arises because the sliding filament model is constrained by more 
than just the inextensibility of the filaments. Its constraints are expressed in the form of 
the strain tensor, or equivalently in the form of the metric tensor in Lagrangian coordinates, 
Eq. (jSJ). It is to be expected that these additional constraints would have consequences for 
the dynamics, just as inextensibility does, through the Lagrange multiplier A. The dynamical 
form that this constraint takes is a certain definite coupling of shear and twist, the last terms 
in Eqs. ( 1591) - (1601) . and more generally every reference to the Lagrange multiplier v and the 
twist velocity u. 

Hines and Blum, in developing detailed nanoscopic models within the framework of the 
sliding filament model, noticed the possibility of such a coupling and wrote more than one 



12 



paper about it 9] 10], pointing out also experimental observations of the effect. It might 
have seemed to their readers, however, that this coupling depended upon nanoscopic details, 
and that there was no necessity for it. No one seems to have pointed out that the coupling 
of shear and twist is a model independent consequence of the sliding filament constraints, 
and must be included in any valid sliding filament model, including dynamical models of 
behavior at the mesoscale. 

APPENDIX A: CILIARY FILAMENTS AS SPACE CURVES 

The rates of change of the vector fields e 2 , e$, along the integral curves of e±, (the filaments 
of the ciliary space) are given by the Lie derivatives (Lie brackets) 

d x e 2 = [ei,e 2 ] = -S x e x (Al) 
d x e 3 = [ei,e 3 ] = -T x e x (A2) 

where the subscript x indicates d x , the derivative along a filament. Thus the curvature 
vector, in the sense of the Frenet equations, is 

kN = S x e 2 + T x e 3 (A3) 

where the curvature is the norm of this vector 



« = v^ 2 + T i ■ ( A4 ) 

The rate of change of the principal normal N along the filament, dotted with the binormal 

B = Sx63 - Tx62 (A5) 

K 

is the torsion 

r = g([e l7 N],B) = T ** S * - S ** T * ( A6 ) 

K 



APPENDIX B: CONSEQUENCES OF FLAT GEOMETRY 

Differential forms dual to the orthonormal vector fields ei, e 2 , e 3 of Eqs. fl9l)- f|TT|) are 

a 1 = dx + Sdy + Tdz (Bl) 
a 2 = dy (B2) 
a 3 = dz. (B3) 
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Define 

A= l -(TS X -ST X + T y -S z ). (B4) 
It is straightforward to verify that the connection forms 

u\ = u\ = cu 3 3 = (B5) 



LU n = — LU 



CO 



x = Sxcr 1 + Aa 3 (B6) 
\ = -lu\ = -T x a l - Aa 2 (B7) 



2 S 
^3 = "Wj 



-Aa 1 (B8) 



satisfy the Cartan structure equations 

da* + Aa 3 =0. (B9) 

Now because the metric is flat, the curvature forms 

9 i j = dJj + J k A u\ (BIO) 

must vanish identically. In particular, 

8\ = -2 A a 2 A a 3 + other components (Bll) 

so A must vanish, and this is Eq. (1271) . The remaining identities reduce to dio\ = duj 1 ^ = 0, 
and this is Eq (|26|) . 

APPENDIX C: GEOMETRICAL CALCULATIONS FOR CILIARY FLOWS 

Under a ciliary flow V it is necessary that gn keep the constant value 1, so that 

= d t g n = 2g{[d x , V], d x ) = 2(a x + Sf3 x + T lx ) , (CI) 

and this is Eq. (|3T|) . To preserve the form of g it is necessary that 

d t g 22 = 2SS t = 2Sd t9l2 (C2) 
d t g 33 = 2TT t = 2Td t g 13 (C3) 
<%23 = 5T 4 + T5 4 = Sd t9l3 + Td t9l2 (C4) 

and these conditions, using Eq. (jlj), are Eqs. (132|) - (!M|) . Finally, it is not enough that g keep 
the form of Eq. (jSJ), it must also continue to have zero Riemannian curvature. The most 
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efficient way to handle this computation is to use the moving orthonormal frame associated 
with the changing g. The change in these vector fields is computed by the Lie derivative 



(the Lie bracket for vector fields) . Thus under the ciliary flow V of Eq. (15U|) 

d t ei = [ei, V) = a x e x + (3 x d y + j x d z = (3 x e 2 + i x e 3 (C5) 

d t e 2 = [e 2 , V] - (d t S) ei = (a\[e 2 , V]) - d t S) e x + (e 27 )e 3 (C6) 

d t e 3 = [e 3 , V] - (d t T) ei = (^([es, V]) - d t T) e x + (e 3 f3)e 2 (C7) 

As a check, one verifies what is required by orthonormality, 

d t S = a\[e 2 ,V)) + p x (C8) 

d t T = a\[e 3 ,V])+ lx , (C9) 

using Eq. (jl]), and also e 3 /3 = — e 2 7 by Eq. flMl) . To summarize, 

9 t e! = (3 x e 2 + 7xe 3 (CIO) 

^e 2 = -p x e x - {e 3 (3)e 3 (Cll) 

d t e 3 = ~lxZ\ ~ (e 2 7)e 2 (C12) 

For the Riemannian curvature to vanish it is necessary that 

= <9i[e 2 ,e 3 ] = [d t e 2 ,e 3 ] + [e 2 ,d t e 3 ] (C13) 

= [(e 3 p) x - {e 2 r r) x + 2T x p x - 2S x ^ x ]e 1 + (e 2 e 2 7)e 2 + (e 3 e 3 /?)e 3 (C14) 

The vanishing of this vector field is conditions Eqs. fl35|) - fl37|) . It remains to show that 
with these conditions the other components of the Riemannian curvature continue to vanish 
identically. Dual to Eqs. flCT0|) - flCT2l) one has 

d t a l = (3 x a 2 + lx a 3 (C15) 

d t a 2 = -^a 1 - (e^a 3 (C16) 

<V 3 = -1.CT 1 - (e 27 )a 2 (C17) 

The Cartan structure equations, which are identities, require 

d t uo\ = (-/3 XX + T^e^a 1 - S x (3 x a 2 - S xlx a 3 (C18) 

d t u\ = {- lxx + S x {e 2l ))a l - T x (3 x a 2 - T xlx a 3 (C19) 

d t u\ = (C20) 
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and a long but straightforward computation then shows that 



fl^dc^. + w'* A w*-) = (C21) 

for all components. 

APPENDIX D: MAKING S AND T CONSTANT IN ONE NORMAL PLANE BY 
REPTATION 



Let a be a reptation, i.e., V = ad x and a x = 0. Then, rearranging Eqs. (!38l)-( !39|) . the 
convective derivative 

S t - aS x = a y (Dl) 

T t - aT x = a z (D2) 

is the rate of change of strain in a fixed normal plane. (The normal planes are not convected 
with the material but are associated with the stationary pattern of the filaments, so to stay 
in a normal plane, subtract the convective term). Thus in a fixed normal plane, since a x = 0, 

(e 2 S) t > = e 2 (S t - aS x ) = e 2 a y = e 2 e 2 a (D3) 
(e 3 T) t > = e 3 (T t - aT x ) = e 3 a z = e 3 e 3 a (D4) 
(e 3 S) t > = (e 2 T) t / = e 2 e 3 a (D5) 

Here the subscript t' means the time derivative at a fixed normal plane, not at fixed x. In 
this normal plane, we can introduce coordinates (y',z') such that e 2 = d/dy', e 3 = d/dz'. 
Then by the Poincare lemma, Eq. (|27|) says that there exists locally a function F(y',z',t) 
such that in that plane S = e 2 F, T = e 3 F. Extend F to be constant along filaments, and 
take a = —F. Then under the reptation V = ad x , 

(e 2 S)tf = -e 2 S (D6) 
{e z T) t , = -e 3 T (D7) 
(e 3 S) f = (ear)* = ~e 3 S = -e 2 T . (D8) 

Thus these quantities decay to zero exponentially with time. In the limit as t — > oo, Eq. ( 1291) 



16 



holds in one normal plane (and then by Eq. (|26|) it holds in all normal planes). 
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